Collective excitation of a Bose-Einstein condensate 
by modulation of the atomic scattering length 
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We excite the lowest-lying quadrupole mode of a Bose-Einstein condensate by modulating the 
atomic scattering length via a Feshbach resonance. Excitation occurs at various modulation fre- 
quencies, and resonances located at the natural quadrupole frequency of the condensate and at 
the first harmonic are observed. We also investigate the amplitude of the excited mode as a func- 
tion of modulation depth. Numerical simulations based on a variational calculation agree with our 
experimental results and provide insight into the observed behavior. 
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Collective excitation of a Bose-Einstein condensate 
(BEC) is an essential diagnostic tool for investigating 
properties of the ultracold quantum state. Fundamen- 
tal information about condensate dynamics can be deter- 
mined from observations of collective modes [l|-[3| , includ- 
ing the effects of temperature , and the dimensional- 
ity of the system 0. In addition, the interplay between 
these modes and external agents, such as random poten- 
tials [i,[9|, lattices [io|,[iii, as well as other atoms 
can be investigated. Examining the excitation spectrum 
of the BEC allows for a detailed comparison with theoret- 
ical models jl7l - [2Qj , an d related quantum systems such as 
superfluid helium jU, [22| and superconductors [23|, [IJ . 

Generically, a collective excitation is generated by the 
modification of the trapping potential of the conden- 
sate [25|. One convenient method is to apply a sudden 
magnetic field gradient, thereby shifting the center of the 
trap and exciting a dipole oscillation about the trap cen- 
ter. One may also suddenly change the curvature of the 
trap to excite quadrupole modes. The lowest-lying m = 
quadrupole mode is characterized by out-of-phase axial 
and radial oscillations. 

If the condensate is not the only occupant of the trap, 
i.e., there exists a thermal component or another species 
of atoms, then the other atoms may also be excited 
through these processes. The evolution of a collective 
excitation can therefore be complicated because the mul- 
tiple components may affect damping or induce frequency 
shifts of the oscillation Therefore modulating the 

trap, although an extremely useful tool for an isolated 
condensate, can be cumbersome when the system to be 
studied is multi-specied. An alternative approach is to 
excite the condensate alone, leaving the other occupants 
of the trap untouched. 

In this work, we demonstrate the excitation of the 
lowest-lying quadrupole mode in a BEC of ''Li by mod- 
ulating the atomic scattering length via a magnetic Fes- 
hbach resonance. In contrast to abruptly changing the 



scattering length [26|, sinusoidal modulation enables the 
controlled excitation of a single mode at a specific fre- 
quency. In addition, by using this method a coexist- 
ing thermal component will be minimally excited by the 
mean- field coupling to the normal gas i28|]. In the 
case of a multi-species experiment, resonant modulation 
of the scattering length of one species will not necessarily 
excite the others, depending on the details of other Fes- 
hbach resonances present in the system [29]. Therefore, 
this technique may be useful for investigating non-zero 
temperature effects, and as a powerful diagnostic tool for 
multi-species ultracold atomic experiments. 

A trapped BEC at zero temperature may be described 
by the three-dimensional cylindrically symmetric Gross- 




FIG. 1: (color online) Quadrupole oscillation excited by 
the modulation of a. A series of in situ polarization phase- 
contrast images of condensates taken during excitation at 
Q = (27r) 10 Hz, separated in time by 15 ms. The change in 
peak density is nearly an order of magnitude from the most 
compressed to the most extended condensate. The horizontal 
scale of these images has been stretched by a factor of two for 
clarity. There is negligible excitation of the dipole mode. 
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Pitaevskii equation [30| 



m 



(1) 



where m is the atomic mass, the trapping potential is 
V = (l/2)m(cj^r^ +cc;^2;^) with uOz (c^r) the axial (radial) 
trapping frequency, a is the s-wave scattering length, and 
the density is given by n = It is convenient to in- 

troduce the anisotropy parameter A = uj^l^r- We use 
a variational approach to solve this equation and deter- 
mine the frequencies of the lowest-lying modes. Using a 
Gaussian ansatz and minimizing the corresponding en- 
ergy functional, we derive the following coupled differ- 
ential equations for the dimensionless axial and radial 
half-widths and Ur of the condensate |mI 



Uz^ \ Uz = ^ 



(2) 



where the interaction parameter P = \/2/7r {Na/lr)^ 
with Ir = ^JhjmbJr the radial harmonic oscillator size, 
and N is the number of condensed atoms. We solve 
Eq. (j2j) in the case of harmonic motion of the Gaussian 
sizes about their equilibrium values. The frequencies of 
the lowest-lying quadrupole oscillation is [3l| 



(1 + A^ - P2,; 



- a/(1-A2 + P2,3)' + 8P, 
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where Pi j = P / (Auq^Uq^) with uqz and uor the equilib- 
rium axial and radial sizes, respectively. The frequency 
of the m = breathing mode is the sum of the two terms 
in Eq. [3] rather than the difference, and is a factor of 
about 60 higher in frequency for our experimental param- 
eters. For the case of a highly elongated trap (A <C 1), 
in the Thomas- Fermi regime (P ^ 1), we find the well 
known relation ujq = ujz ^5/2 [isl, [32| , and in the non- 
interacting limit (P 0), we find ujq 2 ujz-, as ex- 
pected. 

We may also use Eq. (|2]) to determine the time evo- 
lution of the size of the BEG In particular, we are 
interested in the dynamics associated with the modula- 
tion of a using a Feshbach resonance, which has been 
proposed previously [ll, [35|. In a magnetic field P, the 
scattering length near a Feshbach resonance may be de- 
scribed by 



a{B) = asG ( 1 
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where asG is the background scattering length, A is the 
resonance width, and Poo is the resonance location. We 
consider a time-dependent magnetic field of the form 

B{t) = Pav + (^PcOs(m), 



where Q is the modulation frequency. Using this form for 
P, the result of expanding Eq. dU to first order in the 
small quantity SB <C | Pav — Poo | is 



a{t) c:^ ttav + Sacos{Qt), 



(5) 



where 



a{Bg,^) and Sa ■ 



aBG^SB 

(^av - ^oo)^ 



This expression for a is substituted into Eq. ([2]), and a 
fourth-order Runge-Kutta method is used to solve the 
system of equations numerically. The results may be di- 
rectly compared with those from experiment. 

Our methods for creating an ultracold gas of ^Li have 
been described previously [361, • We confine atoms in 
the |1, 1) hyperfine state of ''Li in an optical trap and use 
a set of Helmholtz coils to provide an axially oriented 
bias field. We determine the radial trapping frequency 
by parametric heating to be ujr = (27r) 235(10) Hz, and 
the axial trapping frequency by dipole oscillation to be 
ujz = (27r) 4.85(1) Hz. After evaporation from the op- 
tical trap at 717 G, where a ~ 200 ao, we have TV ~ 
3 X 10^ atoms with a condensate fraction > 90%. We 
then ramp the bias field in 4 s from the initial value to 
Pav = 565 G (where aav ^ 3ao). For these experimen- 
tal values, the dimensionality parameter is AP ^ 0.3, 
close to the transition into the quasi- ID regime (deter- 
mined by AP <C 1) [3^. At this point we oscillate the 
magnetic field with a modulation depth of SB = 14 G, 
corresponding to 2ao. We use in situ polarization 

phase-contrast imaging to obtain the density distribution 
of the condensate [3^] to which we fit a Gaussian charac- 
terized by 1/e axial and radial half- widths. 

In Fig.[TJ we show pictures taken with 1] = (2 tt) 10 Hz. 
A quadrupole oscillation of the cloud is readily observ- 
able. The large oscillation amplitudes considered here ex- 
tend over approximately 1 mm of the optical trap. A har- 
monic approximation of the trapping potential about the 
trap center is less than 10% in error over this range. The 
size of the cloud in Fig. [1] as a function of time is mod- 
eled well by the variational calculation, consistent with 
negligible anharmonic contributions. Furthermore, we 
observe no damping of the quadrupole mode over many 
oscillation periods, consistent with a negligible thermal 
fraction. 

Results from the variational calculation show that dur- 
ing the excitation if ^7 < cjq, then the axial and radial 
sizes of the cloud follow the change in a: growing as a 
increases, and shrinking as a decreases. This in-phase 
behavior of both the axial and radial sizes of the cloud is 
expected for adiabatic changes in a, and therefore should 
not be confused with the high-lying m = breathing 
mode for low frequencies. However when ft > ujq^ the 
radial size follows nearly in-phase, while the axial size 
lags behind the radial size by half a period — an out-of- 
phase oscillation. In both cases when the excitation is 
stopped, the cloud undergoes free quadrupole oscillations 
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FIG. 2: (color online) Axial (red circles) and radial (blue 
squares) 1/e sizes of a condensate during and after excita- 
tion with ttav 3ao and 6a ^ 2ao, where (a) Q = (27r)3Hz 
or (b) Q = (27r) 10 Hz. The natural oscillation frequency is 
ojQ = (27r) 8.2(1) Hz. Note that the 1.8 Hz beat note is the 
cause of the diminished amplitude in (b) after the modula- 
tion drive is turned off. The solid lines are results from the 
variational calculation using the same parameters as in the ex- 
periment and contain no adjustable parameters. The dashed 
lines are representative of the driving excitation. The reso- 
lution of the optical imaging system is ^3 /xm, which limits 
accurate determination of the radial sizes. The frequency of 
the breathing mode is on order 470 Hz and therefore no effects 
of this mode are present. 



with the axial and radial sizes tt out-of-phase. Charac- 
teristic data for these two regimes are shown in Fig. [2] 
and reasonably agree with the variational results. We 
fit the time evolution of the size of the cloud after ex- 
citation and determine the free quadrupolar oscillation 
frequency to be ujq = (27r) 8.2(1) Hz, in good agreement 
with the predicted value of 8.17 Hz from Eq. [3l where we 
have used A = 0.021 and P ^ 15. Similar agreement 
between measured and predicted quadrupole frequencies 
in the dimensional crossover regime have been previously 
observed . The amplitude of this free oscillation is de- 
pendent on the duration of excitation as well as the phase 
of the driving force at the time when the excitation is 
stopped. Therefore, care must be taken when comparing 
data with theoretical predictions of the amplitude during 
the free oscillation. 

A less parameter-dependent measure is to observe the 
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FIG. 3: (color online) Fractional amplitude of the drive (red 
circles) and quadrupole (blue squares) modes during excita- 
tion at frequency Q with a modulation depth ^ 1 ao, where 
aav 3ao. The solid lines are results from the variational 
calculation with no adjustable parameters. The oscillation is 
notably asymmetric for fractional amplitudes of order 1 and 
larger, as shown in Fig. [2jb). 

amplitude of the oscillation during excitation. As can be 
discerned from Fig. O during excitation the size of the 
condensate oscillates at both the drive frequency and the 
natural quadrupole frequency. Beating between these fre- 
quencies modulates the instantaneous deviation from the 
unperturbed size. An ideal method for determining the 
energy in the driven mode and the excited quadrupole 
mode separately is to use Fourier analysis [40|. This 
method can be experimentally difficult, however, given 
the required number of coherent oscillations needed to 
accurately resolve the sinusoidal peak in the Fourier spec- 
trum. Instead, we assume that the system can be de- 
scribed by the linear combination of sinusoids at the 
known frequencies Vt and ujq. After driving an excita- 
tion for 0.5 s, we fit the observed condensate axial size 
to the following expression during an additional 0.5 s of 
excitation: 

u(t) =uo^un sin {Qt^^)^UQ sin {ujqt + (j)) , (6) 

where uq is the equilibrium 1/e size, uq and uq are the 
amplitudes of the drive and quadrupole modes, respec- 
tively, and $ and (j) are the respective phases. The frac- 
tional amplitudes uq/uq and uq/uq are shown as func- 
tions of Q in Fig. [3l As expected, there is a resonant 
enhancement in both the quadrupole and drive modes 
when Q = ujQ. In addition, we see a parametric enhance- 
ment when Q ^ 2ujQ. The dip at exactly 2ujq is due to 
destructive interference between the drive mode and the 
resonantly excited mode. This interference structure is 
observed in both the data and the simulation. 

Larger drive amplitudes push the oscillations into the 
nonlinear regime where the amplitude of oscillation is no 
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FIG. 4: (color online) Quadrupole oscillation driven by a 
large amplitude excitation, 5a ^ Sao where aav ^ 3ao, near 
resonance Q = (27r)9Hz. The non-sinusoidal behavior leads 
to fragmentation of the condensate during the compression 
stages when the axial size becomes smaller than the axial 
harmonic oscillator size of ^17 /j.m. The image at right was 
taken during the compression stage at t — 0.52 s and shows 
fragmentation of the condensate, the field of view is 25 fim x 
380 /xm. 
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FIG. 5: Fractional amplitude of oscillation, when Q ^ cjq, as 
a function of the fractional modulation depth 6a/ aav after 1 s 
of excitation. The solid curve is the result from the variational 
calculation. 



longer linearly dependent on the modulation depth. The 
first noticeable effect is the non-sinusoidal behavior of the 
oscillation seen in Fig.[2fb). In Fig. 2] we show the result 
of driving the system near resonance with a fractional 
modulation depth Sa/a^v = 1. As the amplitude of the 
driven oscillation grows, it eventually becomes compara- 
ble with the original condensate size. At this point, the 
size of the condensate cannot become smaller, and there- 
fore, the oscillation becomes increasingly asymmetric. In 
this manner, the rate of energy transfer into the oscilla- 
tion will decrease. At even larger amplitudes, we observe 
that the condensate appears to fragment [35|. Similar 
looking results have been observed when modulating the 
radial confinement of a cigar-shaped BEC, where Fara- 
day waves may be excited [41|. By fitting the excitation 
spectrum (Fig. [3j) to a skew Lorentzian (where the peak 
excitation frequency, width, amplitude, and skewness are 



fit parameters), we determine the maximum fractional 
amplitude of oscillation of the quadrupole mode after 1 s 
of excitation. Our experimental results for the ampli- 
tude as a function of modulation depth are presented in 
Fig. [5] along with results from the variational calculation, 
which show good agreement with no adjustable parame- 
ters. We note that for the smallest of modulation depths 
investigated, we only observed oscillation of the conden- 
sate when the drive frequency was near resonance. In 
addition, we found a roughly linear scaling of the instan- 
taneous amplitude of the quadrupole oscillation with the 
duration of the excitation in this regime. Furthermore, 
Fig. [5] conveys to us that the data in Fig. [21 which was 
driven at a fractional modulation depth of ^0.7, was not 
in the linear regime, and the data shown in Fig.O driven 
at a fractional modulation depth of ~0.3, had a response 
that departed from hnearity by about 30%. 

Large amplitude oscillations are accompanied by a fre- 
quency shift of the quadrupole mode |42|] . For our geom- 
etry, this shift can be approximated by Suj/ujq ~ O.IA^ 
for small ^4^, where Az = uq/uq is the fractional am- 
plitude of the axial size [40[. The data shown in Fig. [3] 
has Az ^ 2 on resonance, and therefore in this case the 
above approximation for the frequency shift is not valid. 
Our simulated results show a frequency shift ~8% for 
Az = 2 with a shallow amplitude dependence at larger 
Az. Whereas a shift of 10% was observed in the oscil- 
lation of Fig. m we were not able to resolve frequency 
shifts by fitting Lorentzians to the excitation spectra for 
our data at low Az shown in Fig. [5l There is an addi- 
tional frequency shift due to non-zero temperature Q-Q 
which is estimated to be negligible given our low tem- 
peratures and interaction strength [43] . Even though we 
can neglect temperature effects in the measurements pre- 
sented here, our method of excitation of the quadrupole 
mode may be used to study these effects in further detail 
in regimes of stronger interactions. In addition, going 
to large values of a will facilitate investigations of be- 
yond mean-field effects on the collective modes of a Bose 
gas Q , complementary to those observed in a Fermi gas 
at the BCS-BEC crossover EBHiS. 



In this work, we have experimentally demonstrated the 
excitation of the collective low-lying quadrupole mode of 
a dilute Bose gas by modulating the atomic scattering 
length. Our observations are supported by variational 
calculations of the time dependent Gross-Pitaevskii equa- 
tion assuming a Gaussian trial wavefunction. Using this 
formalism we find good agreement with our experimental 
results. Temporal modulation of the scattering length, as 
afforded by Feshbach resonances, provides an additional 
tool for exciting collective modes of an ultracold atomic 
gas. This method is quite attractive in circumstances 
where excitation of the condensate by other means, such 
as trap deformation, is unavailable. In addition, this 
method can be used for condensates in the presence of 
thermal atoms where principal excitation of the conden- 
sate alone is desired, as well as in multi-component gases 
where excitation of only one species can be accomplished. 
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